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Vortices in superconducting strips: interplay between surface effects and 
the pinning landscape 
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Vortices in a narrow superconducting strip with a square array of pinning sites are studied. The interactions of 
vortices with other vortices and with external sources (applied magnetic field and transport current) are calculated 
via a screened Coulomb model. The edge barrier is taken into account and shown to have an important role on 
the system dynamics. Numerical simulations in this approach show that the field dependent magnetic moment 
presents peaks corresponding to history dependent configurations of the vortex lattice. Some effects of the edge 
barrier on the I-V characteristics are also reported. 



1. Introduction 

There is a growing interest in the microscopic 
structure of the vortex lattice (VL) in supercon- 
ducting films with artificially fabricated colum- 
nar defects. In the presence of a perpendicular 
magnetic field and zero current this structure de- 
pends on how strong is the pinning force produced 
by the defects with respect to the internal elas- 
tic forces of the VL. This results in two possible 
situations: commensurate and incommensurate 
lattice, depending on the pinning landscape be- 
ing strong enough to determine the VL structure 
or not. Many experimental jjj and theoretical 
works have reported maxima in the critical 
current j c and field dependent magnetic moment 
m(H) in the points of the phase diagram where 
a commensurate state takes place. 

When a finite current stronger than j c is ap- 
plied, an additional ingredient governs the struc- 
tural evolution of the VL. In the center of mass 
(CM) frame, the vortices interact with a time de- 
pendent pinning potential, which is determined 
by the CM velocity. As the driving current is in- 
creased from j c , a variety of dynamical phases is 
possible. But commensurate ordered lattices are 
only possible at very high driving forces, when the 
pinning potential is averaged in the direction of 
motion becoming a static washboard potential. 
Recently, Carneiro P] has shown by numerical 



simulations that at a high enough temperature 
this ordered state is reached at a finite driving 
force. 

Frequently, these phenomena have been investi- 
gated under the assumption of infinite sample, i.e. 
no surface or geometrical effects are taken into ac- 
count. As it is well known, surface effects can play 
an important and even a dominant role on the ir- 
reversible properties of a superconducting speci- 
men. In the critical state approach, Kuznetsov et 
al. Gj estimated the enhancement due to edge ef- 
fects of hysteresis in magnetization of thin strips 
with strong pinning. The penetration of flux in 
finite sample are also described by the geometric 
barrier model which has shown to be suc- 

cessful in describing the irreversible properties of 
flat high-T c single crystals |7j at high tempera- 
tures. All these continuum models however are 
not able to determine structural details of the VL 
and their relevance in the macroscopic response. 

In the present work, we propose a simple quali- 
tative description of individual vortices dynamics 
in superconducting strips. The vortex-vortex in- 
teractions and the edge barrier are calculated in 
the frame work of the screened Coulomb model 
H . With a change of scale the problem may be 
mapped into the well known problem of an infi- 
nite slab under a parallel field. The calculated 
m(H) presents a series of maxima correspond- 
ing to matching of the VL in the pinning struc- 



2 



ture and a strong assymetry of the upward (vir- 
gin) and downward brenches, which charactere 
the surface effects. We also study the dynam- 
ics of the VL when a transport current is applied 
accross the strip. The strong current inhomogen- 
ity leads to different dynamics in different regions 
of the strip. 

2. Screened Coulomb model and dynamics 

We describe the motion of vortices in a su- 
perconducting strip of thickness d much smaller 
than the London penetration depth A. The ge- 
ometry and coordinate system are depicted in 
Fig. [j]. The vortex equation of motion, ap- 
propriate for overdamped dynamics, is t^v^ = 
Fi(r fc ) - VU p (r k ) + T k (t), where v fc is the in- 
stantaneous terminal velocity of vortex k, Fi(rfc) 
is the total Lorentz Force, U p is the pinning po- 
tential and Tk(t) is a Gaussian Langevin force 
acting on vortex k at instant t. The Lorentz force 
acting on a vortex line is defined by the total 
sheet current J = J jdz probed by this vortex, 
Fl = 0oJ x z. A difficult task is to precisely 
determine an expression for F^ in this system. 
Recently, Brandt has generalized the continuum 
theory to account for the penetration depth and 
numerically calculated the stream lines of current 
of vortices in thin films and washers in the small 
A limit §. 

A usefull procedure is to write the sheet current 
in terms of a scalar potential, the stream function 
g(x,y) satisfying the zero divergence of the cur- 
rent J(x, y) = V x zg = — z x Vg. Note that 
there is a close relation between the stream lines 
probed by a vortex and the energy of this vor- 
tex J! . This is clear when we rewrite the Lorentz 
force as Fl = 0oJ x z = — o Vg. With the help 
of the Biot-Savart law, the London equation for a 
vortex located at rfc in the xy plane of a thin strip 
under a perpendicular field H may be written as 



AX7 2 g(r)-J d 2 r'K{r,r')g(r') = H- 



00 

Mo 
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where K(r,r') is a kernel with a suitable cut- 
off to account for the finite sample thickness ||. 
We propose a simple way of calculating g, us- 
ing a screened Coulomb model approach, suitable 
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Figure 1. As calculated magnetic moment m vs. 
magnetic field H for the thin strip with pinning 
strength Uo = 0.5e. m is in units of (f> and H 
in units of 4>o/[IqA 2 . Inset: system geometry and 
coordinate system. 



for describing thin-film vortices with A compara- 
ble to the system dimension ||. In this large-A 
limit, the screening length A c is identified with 
the smaller of A and some linear dimension of the 
sample. Therefore, for thin strips with W > A, 
we use a nonhomogeneous Helmholtz equation in- 
stead of solving the exact, although much more 
complicated, Eq. [I]. Using the boundary condi- 
tion g = at the edges, this problem is similar 
to that of a superconducting slab under a par- 
allel field 0. Here, we add periodic boundary 
conditions (with periodicity L) in the y direction. 
Solving for a single vortex, one obtains the vortex- 
vortex interaction energy: 



E cosh a m (\yj - y } \ - L/2) 
— — 
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where a m = y I + (rrm/W) 2 , e = 0q//z o A c W. 
The tilded length quantities are normalized by A c . 
With a suitable cutoff to account for the vortex 
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core, one obtains the vortex self-energy: 
cosha m (£ — L/2) . 2 miTXi 
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self 
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(3) 



where £ is the coherence length. Now, solving for 
the external field H, we have 



V H = fo- 



cosh(i 4 - W/2) 
cosh W/2 



- 1 



(4) 



Equations (||) and (Q) form the edge barrier for 
flux entry. Thus, the Lorentz force on a vor- 
tex i due to a vortex j and an external field 
H, for instance, is given by — ViVfo,!^) where 
V(ri,rj) = K^O",,^-) + Keif (ri) + Vff ( r *)- 

3. Periodic pinning: static aspects 

In the above framework, we simulated m(H) 
measurements on a thin strip with W — 3A. 
A square lattice of columnar defects (SLCD) is 
represented by a periodic arrangement of Gaus- 
sian wells with radius r p — 0.045A and period 
a p — 0.3A. The depth of the wells is varied from 
C/o = to (7o = f. The simulations are carried 
out using the same procedure as in Ref. |l(J: at 
each time step a vortex is attempted to nucleate 
in some random point in a line which is a dis- 
tance £ from one of the strip edges. Here we take 
£ = 0.01. The magnetic moment was calculated 
using m = | J d 2 rr x J = | J cPrr-Vg, which for 
the chosen gauge of g = at the edges is simply 
fd 2 rg. 

In Fig. |l| is shown the first m(H) half cycle for 
Uq = 0.5e. This curve depicts features of both 
bulk pinning due to regularly arranged deffects 
and edge effects. Matching peaks are present, 
but behave differently in the upward (virgin) and 
downward branches. As shown in Fig. g, these 
peaks correspond to hystory dependent states of 
commensuration. In upward field, the edge bar- 
rier imposes a vortex free region which produces 
vacant pinning sites even in the second matching 
configuration. In downward field, these states are 
similar to those expected for infinite samples, no 
vortex free region is present. This points up the 
different dynamics of flux entrance and exit. 
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Figure 2. Snapshots of the vortex lattice (closed 
circles) in the matching peaks depicted in Fig. [I]: 
(A) H = 14.5 (upward branch), (B) H = 10.5 
(downward), (C) H = 23 (upward), and (D) H = 
20.5 (downward). The hollow circles represent 
the pinning lattice. 



4. Periodic pinning: dynamical aspects 

We have also simulated measurements of the 
time averaged CM velocity of the VL as a trans- 
port sheet current J t was varied. The current 
is apllied along the y direction and the contacts 
are considered to be far away from the region 
studied, in such a way to guarantee zero current 
divergence in this region. The current distribu- 
tion in the strip width may be calculated from 
the Helmholtz equation for gj with the condi- 
tion that gj must be an odd function and that 
9j(W) —gj(0) — JtW (current continuity). Thus 
the potential energy due to Jt is 



Vj 



■Wsinh^ - W/2) 
2sinhW72 



(5) 



The simulation procedure is the same as de- 
scribed in the preceeding section. We divide the 
strip in three regions of width A c . We computed 
the time-averaged CM velocity vcm in each re- 
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Figure 3. Time-averaged CM velocity vcAi(Jt) in 
the central third of the strip for pinning strengths 
U a = 0,0 .25, 0.5, 0.75, and l.Oe (from top to bot- 
ton). Panels: density-density correlation func- 
tions for Uq = e and different current values (ar- 
rows). 



gion as a function of the transport current after 
the field has been ramped up to 24(/) /^oA^. In 
Fig. [| we show the current dependence of vcm in 
the central region for different pinning strengths 
varying from — to Uq = e. As expected, 
there is a finite critical current for zero pinning, 
since the edge barrier pins the VL. For nonzero 
pinning, two regimes of motion are identified: in 
the low current regime the interstitial vortices are 
depinncd and start moving; at higher value of Jt 
all the lattice is moving. As shown by the density- 
density correlation function of the vortices posi- 
tion (panels of Fig. ||), no reordering occur in 
the current range studied. Reordering would be 
expected at least for weak pinning since the tem- 
perature used in the simulations is a small frac- 
tion of the zero-drive melting temperature. But 
here, lattice ordering is contaminated by the ran- 
dom nature of flux penetration through the sur- 
face barrier. However, long range time order was 
observed. 



5. Summary 

We have described the dynamics of vortices in 
narrrow superconducting strips in the screened 
Coulomb model approach. We calculated the 
vortex-vortex interaction energy and the interac- 
tion with an edge barrier. In this scenario, we 
carried out simulations on vortex motion in field 
ramp measurements of the magnetic moment and 
I-V characteristics using Langevin dynamics. The 
m(H) curves evidenced the strong influence of the 
edge barrier in the VL configuration in upward 
field, where the first matching states are not fully 
accomplished in the entire sample. For the I-V 
characteristics, we observed that the edge barrier 
represents an effective pinning for the VL as ex- 
pected, and ordered motion of the VL is disturbed 
by the random penetration of vortices. A more 
detailed study will be reported elsewhere. 

This work was sponsored by the Brazilian agen- 
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